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Abstract 


We apply analytic perturbation theory to the inclusive decay of a r lepton into 
hadrons. It is shown that the resulting analyticity of the coupling constant 
strongly influences the value of the QCD A-parameter extracted from the 
experimental data on r decay. 
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I. INTRODUCTION 


The inclusive character of the decay of a r lepton into hadrons and the fact that the 
nonperturbative QCD contributions to this process are very small ]l[ make it possible, in 
principle, to describe this process on the basis of the standard methods of quantum field 
theory without any model assumptions. Measurement of the ratio of hadronic to lcptonic 
r decay widths, i.e., the quantity R r = T(r —> hadrons + u)/T(t —> luu) , allows one to 
extract, with a high degree of accuracy, the value of the strong coupling constant «qcd at 
the r mass, Q = M r ~ 1.78 GeV. Comparing this value with the values of «qcd found 
at higher energy, for example, «qcd(Q 2 = M§), is an important test of the applicability 
of QCD perturbation theory over a wide range of energies and requires a careful check. 
At present, both experimental and theoretical investigations of the decay of a r lepton are 
continuing intensively (see the recent reviews in 0). 

One usually employs analytic properties of the hadronic correlation function in order 
to rewrite the original expression for the r hadronic rate, which involves integration over 
a nonperturbative region of small momenta, in the form of a contour integral over a circle 
of sufficiently large radius Q 2 = to apply perturbation theory (PT) p].[Z]. However, the 

perturbative approximation, which introduces a ghost pole, violates the analytic properties 
required to use Cauchy’s theorem in this manner. 

In this paper, we will apply analytic perturbation theory (APT) [[|.[|] in which it is 
possible to maintain the correct analytic behavior. Within this approach the correct analytic 
properties of the running coupling are provided by nonperturbative contributions which 
emerge automatically from dispersion relations.^] The analytic running coupling constant 
obtained in such a way turns out to be remarkably stable for the whole interval of momentum, 
and has a universal infrared limit at Q 2 = 0, independent of the value of the QCD scale 
parameter A. 

For our purpose, it is important that, within this approach, it is possible to give a 
self-consistent definition of the running coupling constant in the Minkowskian (timelike) 
region |J. This fact allows us to obtain two equivalent representations for the QCD correc¬ 
tion to r decay, involving the timelike and the spacelike definitions of the running coupling 
constants, respectively. For our numerical estimations we will use APT at 2-loop order. 


II. WHY IS ANALYTICITY IMPORTANT? 


The initial theoretical expression for R r in the case of massless quarks contains an integral 
over timelike momentum s J3J 



( 1 ) 


1 Some other aspects of the dispersive approach have been discussed in [/]] and in the recent 
papers ||. 
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where the range of integration extends down to small s and cannot be calculated in the 
framework of the standard perturbation theory. The method of calculation of R r based on 
exploiting certain analytic properties of the hadronic correlators II(s) allows one to rewrite 
the expression ([!]) by using the Cauchy theorem in the form of a contour integral in the 
complex s-plane with the contour running counterclockwise around a circle centered on the 
origin of radius M 2 : 


R t 




1 

27ri 


s\=M 2 s 


M3 


M3 


( 2 ) 


where 2 ' 


D(q 2 ) = ~q< 


dU(q 2 
d q 2 


In the representation (||) the contour has a sufficiently large radius, and it is possible, 
in principle, to calculate R r perturbatively. However, the transition to the contour repre¬ 
sentation requires certain analytic properties of the correlator. Namely, the correlator n(s) 
is an analytic function in the complex s-plane with a cut along the positive part of the 
real axis. The parametrization of n by the perturbative running coupling constant violates 
these analytic properties |TlTJ|] . It is clear that the difference in the regions of integration in 
the initial expression ([I]) for R r and the expression (||) obtained after applying the Cauchy 
theorem makes it necessary to parametrize n in Eq. (|l|) and D in Eq. (|2|) with different cou¬ 
pling constants. Indeed, a renormalization-group analysis gives a running coupling constant 
determined in the spacclike (Euclidean) region, while the initial expression (jl]) contains an 
integration over timelike momentum and therefore to calculate Eq. (|2|) requires some pro¬ 
cedure of analytic continuation from spacelikc to timelike momentum. To this end, we will 
determine the effective coupling constant in the spacelike region (f-channel), a eff , and in the 
timelike region, a® ff , (using the definition that a = aQCD/4vr) as (where the overbar signifies 
that the coupling constant has the correct analytic properties) 


D(q 2 ) ~ 1 + dia(q 2 ) + d 2 a 2 (q 2 ) H-= 1 + d 1 a eS (q 2 ) , (3) 

Inffi(s) ~ 1 + ria< j (s) + r 2 a 2 (s ) H-= 1 + ridf^s), (4) 


where the first perturbative coefficients d\ and rq are equal to each other: d\ — r\ — 4. 

The dispersion relations for the Zi-function yield the following relation between the ef¬ 
fective coupling constants: 


d eS (q 2 ) = -q 2 f 
Jo 


ds 


af(*), 


> - q 2 Y 

' s+i£ dz 


-eff/ \ f S 16 dZ q , v 

27T1 Js- ie Z 

Now, separating the QCD contribution A r in R r , 


RT = S(\V ud \ 2 + \V us \ 2 )S EW (1 + A r ), 


( 5 ) 

( 6 ) 

(7) 


2 Here, we will use the standard definition q 2 


—Q 2 , so that in the Euclidean region Q 2 > 0. 
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where V u d and V us are the CKM matrix elements, and Sew is the electroweak factor (see [1]), 
we obtain from Eqs. (1) and (2) the two equivalent representations, 


and 


A r 


2 r 1 


i-m? ds 

o M 2 




dfi 


s) » 


A r 


At dz/ 
2ni J\z\=M% z \ 



( 8 ) 

( 9 ) 


It should be noted that the equivalence of these formulae holds only in the case of the 
above-mentioned analytic properties of the correlator II(s) and the /^-function and that these 
analytic properties are broken in the standard perturbation theory. The QCD contribution 
represented in the form (|2|) is the expression which one usually uses for theoretical analysis. 
In principle, the expression ([)]) can be calculated on the basis of perturbation theory, but 
then it is impossible to return from Eq. (|9j) to Eq. (|8]) which corresponds to the initial 
Eq. ([[]), and nothing can be said about the error associated with switching from Eq. (0) to 
Eq. (|9|). Therefore, it is impossible to give in the framework of the standard perturbation 
theory a self-consistent description of the inclusive decay of a r lepton into hadrons. 


III. THE APT METHOD FOR r-DECAY 


Now, we apply analytic perturbation theory (APT), which was recently proposed 
in PHII, to describe the inclusive decay of a r lepton into hadrons and to give a quantitative 
estimate of the effect associated with Q 2 -analyticity of the running coupling constant. The 
method of APT allows one to avoid the above-mentioned difficulties and evaluate both the 
initial integral over the physical region and the contour representation which are equal due to 
the Cauchy theorem. In the framework of APT, If and the D-function can be parametrized 
by the running coupling constant with the correct analytic properties as in Eqs. @ and ([|), 
where the running coupling constants in the space- and timelike regions can be expressed in 
the terms of the spectral density g(a) as 0 


d{z) = - f 

7T JO 
1 


da 


a - 

00 da 


Q{°), 


i da 

a s [s) = - — g{a ). 

7T J s (J 

In the leading order, the expression for a(z) has the form 


a (1) (A = 


A) 


+ 


ln(— z/A 2 ) 1 + z/A 2 


( 10 ) 

( 11 ) 

( 12 ) 


where /3q = 11 — 2nf/3 is the first coefficient of the f3 function. The first term in (12) 
determines the asymptotic behavior at large momenta and has the same form as in per¬ 
turbation theory. The second term, which appears automatically, reproduces the correct 
analytic properties, and the ghost pole at z = —A 2 does not arise. We underscore once 
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again that APT makes it possible to implement the correct transition from expression (HD to 
expression (H). These expressions simply coincide, as they should, while the application of 
perturbation theory with the standard renormalization-group refinement runs into serious 
difficulties. 

The fundamental quantity in the APT approach is the spectral density g(a) by which 
one can parametrize both the running coupling constants in the spacelike and in the timelike 
regions. We now find a formula which expresses the strong interaction contribution to R T 
via the spectral density function. To this end, introduce an effective spectral density p eS (a) 
that corresponds to the effective coupling constant [see Eq. (]3|)] and write down the effective 
coupling constants in the Euclidean and physical regions as follows, 


d eS (z) 

af(s) 


1 

7T 
1 

7T 


da 


a — 
00 da 

— P 
a 


z 

eff 



(13) 

(14) 


First, consider the contour representation (9j) for the QCD contribution A T . By using 
Eq. (13) this formula can be rewritten as follows 



Of course, the same result ( |T5| ) will be obtained if we substitute Eq. (|T4| ) into Eq. (^). 

Due to the feature of universality, the first term in Eq. (|15|) can be expressed in terms 
of only the first /3-function coefficient. To demonstrate this here let us use the following 
asymptotic representation for the effective coupling constant 


A PT 


(<? 2 ) = 


oo k 

y/ y/ a k,m. 

k =0 m=0 


In" 


t' k 


(16) 


where i = ln(Q 2 /A 2 ) and a 0 ,o = 1/An Calculating the effective spectral density as the 
discontinuity of the effective coupling ([R]), for the infrared limiting value of the analytic 
effective coupling constant we get, with L = hitr/A 2 , 


with 


a AFr(0) 


1 

- dL p eS (L) 

7r J— oo 
-i oo k 

j+EE a t , m Aaf m (0) 

Po fc =1 m= 0 


A 4 ff m(0) 



dL 


In m [L - hr) 
(L - m) k + l ' 


(17) 


(18) 


By repeatedly integrating by parts, it is now very easy to show that 




( 0 ) = j—A. a 

y ’ k m 


—efF 


fc,o(0) — 


0. 


(19) 
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Consider the 2-loop level, in which the perturbative effective running coupling constant 


is 


4t(9 2 ) = a PT (q 2 ) + jf4 T (g 2 ) , 

ttl 


( 20 ) 


where for rif active quarks d 2 = 16(1.9857 — 0.1153n/) [|ll[]. By using Eq. (]20|) we can 
obtain the corresponding effective spectral density. It is important to appreciate that in the 
framework of the APT approach there is very little sensitivity to the approximation used 
in solving the renormalization group equations. This fact is demonstrated in Fig. |Ij. where 
we plot three spectral densities. The solid line corresponds to the spectral density obtained 
from the exact quadrature of the 2-loop (3 function, 


M = f 


a dx 
(3(x) 


1 

a 


A )B\ In 1 + 


PoBia 


( 21 ) 


where B\ = /?i/(/?o) 2 , and A = 102 — 38n//3 is the two-loop coefficient of the /3-function. 

The dashed line corresponds to the spectral density obtained from the asymptotic rep¬ 
resentation of the perturbative 2-loop running coupling, 


a = 




1 ~B 1 


In A 


The dot-dashed line corresponds to the following perturbative running coupling 

1 


d{Q 2 ) = 


( 22 ) 


(23) 


AA + PqB i ln(l + £/ Bi) 

This expression, which we will use for our numerical calculations, is the result of exact 
integration of the two-loop differential RG equation (^1|) solved by one iteration. For this 
case the spectral density which is associated with the running coupling (|23|) is 


Q{v) = 


I(L) 


A R 2 {L) + I\L ) 


(24) 


with 


A 


R(L) = L + B 1 In 
I(L ) = 7r + Bi arccos 


L 

1 + A 
B\ 


+ 
- L 


IT 


(25) 


(Bi + L) + 7T 2 
As a result, for the effective spectral density, we find 

1 d 2 2 R(L)I(L) 


P° ( ff ) = q(o) + 


PUA&m+pa)}- 


(26) 


By substituting this formula into Eq. (|15D and using the fact of universality for the first 
term we obtain the QCD correction for the inclusive decay rate of the r lepton. For our 
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0.03 



L = ln(a/A 2 ) 

FIG. 1. Plot of the spectral densities. The solid line refers to the exact quadrature of the RG 
equation, the dashed line refers to the asymptotic representation from the perturbative running 
coupling, and the dot-dashed line corresponds to solution obtained by one iteration. 


numerical estimations we will use the world average R T = 3.633 ± 0.031 [O, which leads to 


the following values of the APT running couplings: a s (M%) = 0.378 ± 0.026 in the timelike 
region and a(M^) = 0.400 ± 0.026 in the spacelike region, respectively. The reader should 
carefully note that the values of the strong coupling constant in the timelike and the spacelike 
region do not agree, even though they are evaluated at the same scale. This difference could 
become significant experimentally as the precision of the measurements improves. 

It should be stressed that these values of the coupling constants have been obtained 
with a value of the QCD scale parameter much larger]] than that found in conventional 
perturbation theory: for three active flavors, A( 3 ) = 935 ± 159 MeV. We have found that 
the value of A is very sensitive to the value of R r . For example, if we were to use a smaller 
value of R r = 3.559 ± 0.035 ||14|| , we would obtain A = 640 ± 127 MeV. Note here that, 


3 The same observation has been made in the one-loop level in [13|. 
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as we have demonstrated, the conventional perturbative parameterization of the D-function 
is inconsistent with the required analyticity properties and it is not possible to write the 
contour integral (|2|) in order to extract Aq CD . Nevertheless, in order to give the reader a 
feeling for the relation between the APT and PT values of A, we require that that the 2-loop 
D-functions be the same in both schemes at the r scale, and obtain the relation between 
the two definitions of A shown in Fig. |2j. 



FIG. 2. Aapt vs. Apt at the r-mass scale. 


It is possible to understand in a simple manner why in this approach the A parameter 
is so much larger than the PT value. To this end, we will use the approximate formula for 
the 2-loop APT coupling constant 


a(q 2 ) 




a PT(q 2 ) + 


l l 

W 1 + q 2 / A 2 


Cl A 2 

+ Po q 2 


(27) 


where, for three active quarks, C\ = 0.035. The accuracy of this formula on the interval 
2.5 < |g|/A < 3.5 is about 0.5%, and Eq. (^7j) practically coincides with the exact formula 
for larger values of momentum. A key feature of this formula is that the perturbative and 
nonperturbative contributions are separated, as in the 1-loop case, which allows us to write 
the QCD correction to r decay as the sum of perturbative and nonperturbative terms. The 
main part of the corresponding nonperturbative contribution is 
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ck 

A) 


(28) 


Anp 


1 + 2 


.d2 V 
di f3 0/ 


M? ' 


Due to the negative sign of this contribution, the perturbative term must be larger than in 
the case of PT (which implies a larger value of A), in order to obtain the same value of the 
QCD correction. 


IV. CONCLUSION 

We have considered the method of analytic perturbation theory and its application to 
the semilcptonic decay of the r lepton. It should be stressed that the principle of causality, 
the consequence of which are certain analytic properties of the running coupling constant, 
leads to the essential stability of the running coupling constant in the infrared region with 
respect to higher loop corrections. Here, the prime point is the universal value of the analytic 
coupling constant at Q 2 = 0, which does not depend on the experimental estimates of the 
QCD scale parameter, A, nor on the number of loops in which we approximate the running 
coupling constant. In other words, families of curves corresponding to different values of 
A and to different numbers of loops have the common point a(Q 2 = 0) = 47 t//3 0 , and are 
represented by a bundle of curves. 

Moreover, this approach allows us to give a self-consistent definition of the running 
coupling constant in the timelike region. Thus, in this approach it is possible to give two 
equivalent expressions for the QCD correction to inclusive r decay, either in terms of the 
timelike coupling, or in terms of the coupling defined for complex Euclidean momentum. 
This is not possible in conventional perturbation theory, due to the breaking of analyticity. 

The fact that we are now able to define a consistent timelike coupling constant enables 
us to perform matching between regions with different numbers of active flavors in the phys¬ 
ical region, where at least in the leading order the number of active quarks is an obvious 
fact associated with physical thresholds. As a result of this matching procedure the cou¬ 
pling constant in the Euclidean region will know about all physical thresholds through the 
dispersion relation. 

Analytic perturbation theory, in effect, incorporates power corrections (nonperturbative 
effects) into perturbation theory in order to secure the required analytic structure. We may 
refer to these as short-distance perturbative power corrections 0. There are indications 
that the ambiguities connected to the asymptotic character of the conventional perturbative 
expansion and the ambiguities in the nonperturbative matrix elements should be similar to 
each other. The convergence properties of the new expansion are different from those of the 
conventional expansion, and empirically the APT seems to be convergent. Therefore, the 
role of the nonperturbative power corrections, which describe the long-distance dynamics 
within the standard approach (in the language of the operator-product expansion, they are 
associated with quark and gluon condensates), is undoubtedly changed. We will consider 
the importance of long-distance nonperturbative power corrections in this scheme elsewhere. 

The results of the analysis performed above demonstrate the importance of analytic¬ 
ity in the running coupling constant, not only from the fundamental point of view—a 
self-consistent theoretical description of r decay—but also from the standpoint of giving 
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a self-consistent description of the Q 2 evolution of the coupling constant and extracting the 
parameter Aq CD from the experimental data on r decay. 
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